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Abstract 

In a recent publication we have calculated elastic charmonium production in ep 
collisions employing realistic charmonia wave functions and dipole cross sections and 
have found good agreement with the data in a wide range of s and Q 2 . Using the 
ingredients from those calculations we calculate exclusive electroproduction of charmo- 
nia off nuclei. Here new effects become important, (i) color filtering of the cc pair on 
its trajectory through nuclear matter, (ii) dependence on the finite lifetime of the cc 
fluctuation (coherence length) and (iii) gluon shadowing in a nucleus compared to the 
one in a nucleon. Total coherent and incoherent cross sections for C, Cu and Pb as 
functions of s and Q 2 are presented together with some differential cross sections. The 
results can be tested with future electron-nucleus colliders or in peripheral collisions of 
ultrarelativistic heavy ions. 
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1 Introduction 



In contrast to hadro-production of charmonia, where the elementary production vertex is still 
debated, electro(photo)production of charmonia seems well understood: the cc fluctuation 
of the incoming real or virtual photon interacts with the target (proton or nucleus) via the 
dipole cross section a qq and the result is projected on the wave function of the observed 
hadron. In the 7*p — > reaction only the dipole cross section a qq and the wave function 
of the \I> enter stands for J ftp or ip f ). In comparing calculations with experiment (which 
fortunately are rather precise and cover a wide range of values s and Q 2 of the incoming 
7*), various parameterizations of the dipole cross section and the wave function have been 
successfully tested |J. In the exclusive electroproduction of charmonia 7* A — > ^/X, where 
X = A (coherent) or X = A* (incoherent, where A* is an exited state of the A-nucleon 
system) new phenomena are to be expected: color filtering, i.e. inelastic interactions of the 
cc pair on its way through the nucleus is expected to lead to a suppression of \I> production 
relative to Ax 7 . p _>^ p . Since the dipole cross section a qq also depends on the gluon distribution 
in the target {p of A), nuclear shadowing of the gluon distribution is expected to reduce a qq 
in a nuclear reaction relative to the one on the proton. Production of a cc pair in a nucleus 
and its absorption are also determined by the values of the coherence length l c and the 
formation length If. 

Explicit calculations have already been performed in in the approximation of a short 
coherence (or production) length, when one can treat the creation of the colorless cc pair as 
instantaneous, 

where v is the energy of the virtual photon in the rest frame of the nucleus. At the same 
time, the formation length may be long, comparable with the nuclear radius Ra, 

2v 

l f = M 2 _ M 2 ~ r a ■ (2) 

In the wave function formation is described by means of the light-cone Green function 
approach summing up all possible paths of the cc in the nucleus. The result has been 
unexpected. Contrary to naive expectation, based on the larger size of the ip' compared to 
J/ifj, it has been found that if)' is not more strongly absorbed than the Jfifi, but may even be 
enhanced by the nuclear medium. This is interpreted as an effect of filtering which is easy to 
understand in the limit of long coherence length, If ^> Ra, when the ratio of cross sections 
on nuclear and nucleon targets take the simple form. Indeed, the production rate of ip' on a 
proton target is small due to strong cancellations in the projection of the produced cc wave 
packet onto the radial wave function of the ip' which has a node. After propagation through 
nuclear matter the transverse size of a cc wave packet is squeezed by absorption and the 
projection of the ip' wave function is enhanced |||| since the effect of the node is reduced. 

However, the quantitative predictions of are not trustable since the calculations have 
been oversimplified and quite some progress has been made on the form of the dipole cross 
section a qq and the light cone wave functions for the charmonia. Therefore we take the 
problem up again and provide more realistic calculations for nuclear effects in exclusive 
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electroproduction of charmonia off nuclei relying on the successful parameter free calculations 
which have been performed recently by the authors ]I| for elastic virtual photoproduction of 
charmonia, 7* p — ■> \& p. Since the present paper is a direct follow up of [|lj we do not repeat 
all definitions and notations, but rather request the reader to consult [0] for more details. 

Whenever one deals with high-energy reactions on nuclei, one cannot avoid another prob- 
lem of great importance: gluon shadowing. At small values of x, gluon clouds overlap in 
longitudinal direction and may fuse. As a result, the gluon density per one nucleon in a 
nucleus is expected to be reduced compared to a free proton. Parton shadowing, which 
leads to an additional nuclear suppression in various hard reactions (DIS, DY, heavy flavor, 
high-j>r hadrons, etc.) may be especially strong for exclusive vector meson production like 
charmonium production which needs at least two gluon exchange. Unfortunately, we have 
no experimental information for gluon shadowing in nuclei so far, and we have to rely on the 
available theoretical estimates PH3]. 

In the present paper we work in the approximation of long coherence time l c 3> Ra 
since calculations can be done with great confidence. Only recently |9j the light-cone (LC) 
dipole approach has been generalized also for the case of a finite coherence length, relevant 
for the vast majority of available data. However, because of technical difficulties those 
calculations employ the oc r\ shape of the dipole cross section. It is still a challenge to 
perform calculations with a realistic dipole cross section for the regime of finite coherence 
length. Nevertheless, we make corrections for finite values of l c employing the approximation 
suggested in [|J. 

The predictions in this paper are meant as a realistic basis for the planning of future 
experiments for electron-nucleus collisions at high energies like in the eRHIC project. It 
would be highly desirable to have good data in order to test the assumptions and ingredi- 
ents which have entered the calculation and which are important aspects of modern strong 
interaction many-body physics. Furthermore, progress in the understanding of charmonium 
electroproduction will also have important consequences for hadroproduction. 

The paper is organized as follows. Section § introduces some notations and the most 
important definitions and expressions. Then the integrated cross sections for coherent and 
incoherent charmonia production are calculated as function of s and Q 2 of the virtual photon 
and the results are presented as a ratio to A-times the cross sections on a proton. In the 
following two sections, two important modifications are introduced. The effect of gluon 
shadowing in the dipole cross section (Section |3|) and the case of finite values of l c {l c ~ 
Ra) (Section f|). Appendix A describes some very useful calculational procedures, which 
considerably increase the accuracy and the speed of numerical calculations. 

2 Shadowing and absorption for cc pairs in nuclei 

Exclusive charmonium production off nuclei, 7*A — > is called coherent, when the nucleus 
remains intact, i.e. X = A, or incoherent, when X is an excited nuclear state which contains 
nucleons and nuclear fragments but no other hadrons. The cross sections depend on the 
polarization e of the virtual photon (in all figures below we will imply e = 1), 

a^ A (s, Q 2 ) = a^ A ( s , Q 2 ) + e a^ A (s, Q 2 ) , (3) 
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where the indexes T, L correspond to transversely or longitudinally polarized photons, re- 
spectively 

The cross section for exclusive production of charmonia off a nucleon target integrated 



over momentum transfer pi is given by 



(7; 



S ,Q 2 



* K 9 -(r T ,s)|7, 



T,L 

cc 



(4) 



where \l/(fY,a) is the charmonium LC wave function which depends on the transverse cc 
separation ry and on the relative sharing a of longitudinal momentum [1]]. Both variables 
are involved in the integration in the matrix element Eq. (|]). \l/(rr, a) is obtained by means 
of a Lorentz boost applied the solutions of the Schrodinger equation. This procedure involves 
the Melosh spin rotation Jll],|l^] which produces sizable effects. 

In Eq. (§ 7 c g' L (f T , a, Q 2 ) is the LC wave function of the cc Fock component of the photon. 
It depends on the photon virtuality Q 2 . One can find the details in Ref. [p]] including the 
effects of a nonperturbative qq interaction. 

The cross section a qq (rT, s) describes the interaction of a colorless quark- ant iquark dipole 
of separation tt and cm. energy squared s with a nucleon. At small values of vt the energy 
dependence should come in the combination x ~ l/(r^s) in order to respect Bjorken scaling. 
A simple parameterization of a qq (r T ,x) suggested in (T^j (GBW) well describes the data 
for F2(x,Q 2 ) at small x and high Q 2 . At low Q 2 , however, the energy variable s is more 
appropriate. A parameterization for a qq (rT, s) suggested in || (KST) describes DIS data 
only up to Q 2 ~ 10 GeV 2 which is a scale relevant for charmonia. Since we cannot see 
why one parameterization should be better than the other, we calculate results for both 
parameterizations as has been done in [|TJ. 

The cross sections for coherent and incoherent production on nuclei will be derived under 
various conditions imposed by the coherence length Eq. (H). At high energies the coherence 
length Eq. ([[]) may substantially exceed the nuclear radius. In this case the transverse size 
of the cc wave packet is "frozen" by Lorentz time dilation, i.e. it does not fluctuate during 
propagation through the nucleus, and the expressions for the cross sections, incoherent (inc) 
or coherent (coh), are particularly simple 0, 



a 



coh 



[s,Q 2 



d 2 bT A (b) 



a q g(r T , s) exp 



d 2 b 



1 



exp 



--a qq {r T ,s)T A {b) 



■^^qq( r T,s)T A (b) 



7, 



T,L 



7, 



T,L 



(5) 
(6) 



Here T A (b) = /f^ dz p A (b, z) is the nuclear thickness function given by the integral of the 
nuclear density along the trajectory at a given impact parameter b. 

In Eqs. d^) and © we also include a small correction due to the real part of the 7*iV — > 
amplitude via replacement 



qq 



1 - 



in d In a, 



2 .9 Ins 



(7) 



As it is common practice we express nuclear cross sections for incoherent and coherent 
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production in the form of the ratio, 

R*(s,Q 2 ) 



a 



7 *A, 



S ,Q 2 



AaT N (s,Q* 



where the numerator stands for the expressions Eqs. (||) and (|6|) for incoherent and coherent 
cross sections, respectively. The details of the technically rather complicated numerical 
calculations of the integrals in Eqs. (|5|) and (|6|) can be found in |Appendix A| . 

The ratio Eq. @ for incoherent electroproduction of Jfip and ip' is shown in Fig. [I] as a 
function of yfs. We use the GBW JT3] and KST |J parameterizations for the dipole cross 
section and show the results by solid and dashed curves, respectively. Differences are at most 
10 - 20 ' 



70. 



rytnc 




TDtnc 




GeV -s/s GeV 

Figure 1: Ratios R^ c for J/ip and ip' incoherent production on carbon, copper 
and lead as function of y/s and at Q 2 = 0. The solid curves refer to the GBW 
parameterization of a q q and dashed one refer to the KST parameterization. 

Analyzing the results shown in Fig. [I], we observe that nuclear suppression of J/ip produc- 
tion becomes stronger with energy. This is an obvious consequence of the energy dependence 
of <Jqq(rT, s), which rises with energy (see in Ref. |rj). For ip' the suppression is rather similar 
as for J/0. In particular we do not see any considerable nuclear enhancement of ip' which 
has been found earlier 0, lij, where the oversimplified form of the dipole cross section, 
°"gg( r T) oc r\ and the oscillator form of the wave function had been used. Such a form of the 
cross section enhances the compensation between large and small distances in the wave func- 
tion of ip' in the process 7*p — > ip'p. Therefore, the color filtering effect which emphasizes the 
small distance part of the wave function leads to a strong enhancement of the %p' production 
rate. This is why using the more realistic r^-dependence of <J q q{rT) leveling off at large tt 
leads to a weaker enhancement of the ip'. This effect becomes even more pronounced at 
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higher energies since the dipole cross section saturates starting at a value ~ r (s) where 
?"o(s) decreases with energy. This observation probably explains why the if)' is less enhanced 
at higher energies as one can see from Fig. [I]. 

Note that the "frozen" approximation is valid only for l c ^> Ra and can be used only at 
y/s > 20 — 30 GeV. Therefore, the low-energy part of the curves depicted in Fig. [I] should 
be corrected for the effects related to the finiteness of l c . This will be done in Section |4]. 

One can change the effect of color filtering in nuclei in a controlled way by increasing 
the photon virtuality Q 2 thereby squeezing the transverse size of the cc fluctuation in the 
photon. For a narrower cc pair the cancellation which is caused by the node in the radial 
wave function of if)' should be less effective. One expects that the if>' to J/if> ratio on a proton 
target increases with Q 2 , as is observed both in experiment and calculation (Fig. 9 of [|TJ). 
Fig. H shows the result for nuclear targets: for large values of Q 2 ratios Rjfi and Rfj) c become 
very similar. This effect has also been observed in the E665 experiment [15] for exclusive 
production of p and p' off nuclei at high energies. 
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Figure 2: Ratios R™^ and Rfj) c for incoherent production on three nuclei dis- 
played function of Q 2 + Mf at fixed s = 4000 GeV 2 . The lines are for the 
GBW (solid curves) and KST (dashed curves) parameterizations. 

Nuclear effects for ip' production shown in Fig. [| demonstrate an even more peculiar 
Q 2 dependence. The overlap of the produced cc state and the ip' wave function rises with 
Q 2 both in the numerator and denominator of the ratio Eq. (||) due to the node effect. 
However, the nuclear filtering for large size cc pairs is especially strong at small Q 2 , hence 
the Q 2 squeezing does not affect the numerator as much as the denominator. This is why 
the ratio R™ C (Q 2 ) can be a falling function at small Q 2 . At high Q 2 , however, the size of 
the cc is so small that color filtering is not an important effect any more, while the nucleus 
becomes more transparent eventually causing R™, C (Q 2 ) to rise. 
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Cross sections for coherent production of charmonia on nuclei are calculated analogously 
using Eq. (||). The results for the energy dependence are depicted in Fig. |3] and for the Q 2 
dependence in Fig. |]. It is not a surprise that the ratios exceed one. In the absence of cc 



jDcoh 

n .J/f 




10 100 10 100 

y/s GeV GeV 
Figure 3: The ratios Rffy and R^i h for coherent production on nuclei as a 
function of y/s. The meaning of the different lines is the same as in Fig. [I]. 

attenuation the forward coherent production would be proportional to A 2 , while integrated 
over momentum transfer, the one depicted in Fig. |3|, behaves as A 4//3 . It is a result of our 
definition Eq. (|) that R$ h exceeds one. To avoid this artifact one may redefine the ratio 
for coherent processes as, 



K h (s,Q 2 



■y*A, 
°coh • 



s,Q 2 



A (T A ) 16tt Brr**(ff N (s, Q 2 ) (T A ) 16ttB^ n 
where B 1 * 1 ^ is the slope of the 7*iV — > ^>N differential cross section and 

(TA) = jjd 2 bT 2 (b) 



(9) 



(10) 



is the mean nuclear thickness. This ratio approaches one when nuclear effects disappear. 
Nevertheless, we will follow the usual definition Eq. (|8|) which is widely used for the presen- 
tation of data. 

We can also predict the dependence on the momentum transfer kx for the charmonium 
electroproduction on nuclei. In the case of incoherent production this dependence is the same 
as for production on free nucleons. However, in coherent production the nuclear formfactor 
comes into play and one has 
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Figure 4: The ratios R c f^ and R c ^, for coherent production on various nuclei as 



a function of Q . The meaning of the lines is the same as in Fig. 

In Fig. |5| we plot the ratios of total distributions (sum of T and L components Eq. (|TI|) in 
the form of Eq. (0) normalized at Q 2 = and kx = 0), 



n(s,Q 2 ,k T ) 



d^(s,Q 2 ) 



da^(s,Q 2 = 0) 



d 2 ki 



d 2 ki 



:i2) 



k T =0 



as functions of k T at s = 4000 GeV 2 for different virtualities Q of the photon. We see that 
the kx dependences are rather similar for J/0 and ip' . The shape of the distribution is 
governed mainly by the nuclear geometry (and not by the size of the (small) charmonium). 
The calculated curves show the familiar diffraction pattern known from elastic scattering on 
nuclei. 



3 Gluon shadowing 

The gluon density in nuclei at small Bjorken x is expected to be suppressed compared to 
a free nucleon due to interferences. This phenomenon called gluon shadowing renormalizes 
the dipole cross section, 

a q q(r T , x) a q q(r T , x) R G (x, Q 2 ,b) . (13) 

where the factor Rc(x,Q 2 ,b) is the ratio of the gluon density at x and Q 2 in a nucleon 
of a nucleus to the gluon density in a free nucleon. No data are available so far which 
could provide direct information about gluon shadowing. Currently it can be evaluated only 
theoretically. In what follows we employ the technique developed in 0. 
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Figure 5: Ratios TZj/^p and TZ^i as functions of kx at s = 4000 GeV 2 for different 
values of Q. All curves are calculated with the GBW parameterization of the 
dipole cross section a q q. 

The interpretation of the phenomenon of gluon shadowing depends very much on the 
reference frame. It looks like glue-glue fusion in the infinite momentum frame of the nucleus: 
although the nucleus is Lorentz contracted, the bound nucleons are still well separated since 
they contract too. However, the gluon clouds of the nucleons are contracted less since they 
have a smaller momentum fraction ~ x. Therefore, they do overlap and interact at small x, 
and gluons originating from different nucleons can fuse leading to a reduction of the gluon 
density. 

Although observables must be Lorentz invariant, the space-time interpretation of shad- 
owing looks very different in the rest frame of the nucleus. Here it comes as a result of 
eikonalization of higher Fock components of the incident particles. Indeed, the nuclear effect 
included via eikonalization into Eqs. corresponds to the lowest cc Fock component 

of the photon. These expressions do not include any correction for gluon shadowing, but 
rather correspond to shadowing of sea quarks in nuclei, analogous to what is measured in 
deep- inelastic scattering. Although the phenomeno logical dipole cross section cr qq -(x, Q 2 ) in- 
cludes all possible effects of gluon radiation, the eikonal expressions Eqs. fl5|)-(§jD assume that 
none of the radiated gluons takes part in multiple interaction in the nucleus. The leading 
order correction corresponding to gluon shadowing comes from eikonalization of the next 
Fock component which contains the cc pair plus a gluon. One can trace on Feynman graphs 
that this is exactly the same mechanism of gluon shadowing as glue-glue fusion in a different 
reference frame. 

Note that Eqs. (§])-(§]) assume that for the coherence length l c 3> Ra- Even if this 
condition is satisfied for a cc fluctuation, it can be broken for the ccG component which 
is heavier. Indeed, it was found in |M that the coherence length for gluon shadowing as 
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about an order of magnitude shorter than the one for shadowing of sea quarks. Therefore, 
one should not rely on the long coherence length approximation used in (^)-(^), but take 
into account the finiteness of If . This can be done by using the light-cone Green function 
approach developed in |J. 

The factor Rg(x, Q 2 , b) has the form, 



(14) 

where afa*N) is the part of the total 7*iV cross section related to a cc fluctuation in the 
photon, 

1 

afa*N) = J d 2 r T J da |*£(r T , a, Q 2 )\ 2 a m {r T ,x) . (15) 
o 

Here ^>2c far, oe, Q 2 ) is the light-cone wave function of the cc pair with transverse separation 
tt and relative sharing of the longitudinal momentum a and 1 — a (see details in 0). The 
numerator Act (7* A) in flnj ) reads ||, 



A / ik A \ f ,9, /\,,9 d? <T fa* N ~ > CCGN) 



(16) 

<?T=0 



00 00 

X 

—00 —00 



J dzi J dz 2 Q(z 2 - zi) p A (b, z x ) p A (b, z 2 ) exp [—i q L (z 2 - Zi) 



Here the invariant mass squared of the ccG system reads, 

M 2 = yjA±R i (17) 

i a i 

where the sum is taken over partons (ccG) having mass mi, transverse momentum hi and frac- 
tion a,i of the full momentum. The ccG system is produced diffractively as an intermediate 
state in a double interaction in the nucleus. z\ and z 2 are the longitudinal coordinates of the 
nucleons ~N\ and N 2l respectively, participating in the diffractive transition 7* N\ — > ccG N± 
and back ccG N 2 — > 7* A^ 2 . The value of Act is controlled by the longitudinal momentum 
transfer 

Q 2 + M 2 

QL = — ^— , (18) 

which is related to the gluonic coherence length If = 1/qi- 

The propagation and interaction of the ccG system in the nuclear medium between the 
points z\ and z 2 is described by the Green function G cE G{r 2 , p 2 , z 2 ; fi, piZy) , where r*i )2 and 
p 1)2 are the transverse separations between the c and c and between the cc pair and gluon 
at the point z\ and destination z 2 respectively. Then the Fourier transform of the diffractive 
cross section in (|IBD 



in 



d 2 afa*N -»• XN) 
dM 2 v ' 



dM\ dq\ 



cos [q L (z 2 - zx)) (19) 

<?T=0 
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can be represented in the form, 



-Re J d\ 2 d 2 p 2 d\id 2 pi J da q d\n(a G ) (20) 
x F^ ca3 (r 2 ,p 2 ,a q ,a G ) G c5G (f 2 ,p 2 ,z 2 ;r 1 ,p 1 ,z 1 ) F T ^ cSG (r 1 , pi, a q , a G ) . 

Assuming that the momentum fraction taken by the gluon is small, a<l, and neglecting 
the cc separation rCpwe arrive at a factorized form of the three-body Green function, 

GccG{r2,P2,z 2 ;r 1 ,p 1 ,z 1 ) => G C c{f 2 ,z 2 ;f 1 ,z 1 ) G GG (p 2 , z 2 ; pi, z x ) , (21) 

where G GG (p 2 , z 2 ; pi, z x ) describes propagation of the GG dipole (in fact the color-octet 
cc and gluon) in the nuclear medium. This Green function satisfies the two dimensional 
Schrodinger equation which includes the glue-glue nonperturbative interaction via the light- 
cone potential V(p, z), as well as interaction with the nuclear medium. 



i — G GG (p 2 , z 2 ; px, Zx 
dz 2 



2va G (l - a G ) 



G GG (p 2 , z 2 ; pi, z x ) , (22) 



where 

2lmV(p,z) = -a GG (p) p A {b, z) , (23) 
and the glue-glue dipole cross section is related to the qq one by the relation, 

9 

cr GG (r, x) = - a q q(r, x) . (24) 
Following |J we assume that the real part of the potential has a form 

ReV(p,z) = - b jf r . (25) 

2va G (l - a G ) 

The parameter bo = 0.65 GeV was fixed by the data on diffractive gluon radiation (the triple- 
Pomeron contribution in terms of Regge approach) which is an essential part of Gribov's 
inelastic shadowing []T7[]. The well known smallness of such a diffractive cross section explains 
why b is so large, leading to a rather weak gluon shadowing. In other words, this strong 
interaction squeezes the glue-glue wave packet resulting in small nuclear attenuation due to 
color transparency. 

In order to get an analytical solution of Eq. (^3|) we assume that <J GG (r, s) ~ C GG (s) r 2 , 
with C GG (s) = da GG (r, s)/dr 2 \ r=0 . We also use the approximation of constant nuclear 
density p^( r ) = Po@(-Ra ~ r ) which is rather accurate for heavy nuclei. Correspondingly 



T A (b) = p L{b), where L(b) = 2^jR\ - b 2 . 

With these simplifications we can perform integration over z 2 in (|16|) and the gluon 
shadowing correction takes the form, 

d 2 b m 

1 - R G (x, Q 2 ) = plj — J dz (L(b) - z) W(x, Q 2 , b, z) , (26) 
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where W(x,Q 2 ,b, z) reads ||, 

4a s (Q 2 ) 



W(x,Q 2 ,b,z) 



3nC qq (x) 



clolq — -— — ~; Re < exp 
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~(Q 2 + 0.25 GeV 2 ) /A 5 
min(x/a G ,a; max ) , 
Q 2 + M 2 , 
bl + a G Q 2 , 

\jV - ia G (l-a G ) v p C GG (x,b) 

%Bj (a G (l-a G )v) , 

B/b 2 , 

t cosh (Qz) + sinh (Qz) , 

(1 + t 2 ) sinh (Qz) + 2t cosh (Qz) . 



(27) 
(28) 

(29) 
(30) 
(31) 

(32) 
(33) 
(34) 
(35) 
(36) 



Here x max = 0.1 and A = 200 MeV. The factor C GG (x, Q 2 , b) is determined by the condition 
of equivalent descriptions in the limit if 3> Ra using the realistic parameterization for 
^gg^t) and its simplified form a GG (r T ) = C GG r\, 



C GG (x,Q 2 ,b)r 2 T T A (b) 



cr qq (x, r T ) R G (x, Q 2 , b) T A (b) 



Here 



(firr)} 



Kf T |VWIi(i-^ 



Sd?r T \* qG (r T )\lJ(r T ) ' 
and the nonperturbative light-cone wave function of the GG (in fact cc and G) dipole. 



|* GG (r T )| 2 ff oc 



-b 2 r 2 
g T 



(37) 



(38) 



(39) 



Obtaining the gluon shadowing factor R G (b) one can relate it to the corresponding value 
T A (b). This gives R G as a function of thickness T A , R G (T A ). This function is then used in 
calculations where T A (b) is obtained with the realistic Woods-Saxon nucleus density. 

Fig. |6| shows the ratios of cross sections calculated with and without gluon shadowing 



S G (s,Q' 



G 1 



a^ A (s,Q< 



(40) 



for incoherent and coherent exclusive charmonium electroproduction. We see that the onset 
of gluon shadowing happens at a cm. energy of few tens GeV. This is controlled by the 
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same. All curves are calculated with the GBW parameterization of the dipole 
cross section a n 
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nuclear formfactor which depends on the longitudinal momentum transfer q c = l/l^. It was 



found in [16] that the coherence length for gluon shadowing is rather short, 

1 



G 



10 xrriN 



(41) 



where x in our case should be an effective one, x = (Q 2 + M£)/2v. The onset of shadowing 



should be expected at q 2 ~ 3/(R c J/) 2 corresponding to 



s G ~ 10m N Rf(Q 2 + M 2 )/V3 



(42) 



where (R'a) 2 is the mean square of the nuclear charge radius. This estimate is in a good 
agreement with Fig. ^. Remarkably, the onset of shadowing is delayed with rising nuclear 
radius and Q 2 . 

Fig. |7| shows the effect of shadowing for the kx dependence of the coherent photoproduc- 
tion. As one can see these corrections are quite small for Q 2 = (and even smaller with 
increasing Q 2 , see Eqs. and 



4 Finite coherence length. 

A strictly quantum-mechanical treatment of a fluctuating qq pair propagating through an 
absorptive medium based on the LC Green function approach has been suggested recently 
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Figure 7: Ratio TZ^ function of k T at s = 4000 GeV 2 and Q 2 = with 
gluon shadowing (solid curves) and without (dashed curves) with the GBW pa- 
rameterization of the dipole cross section. 

in 0. However, an analytical solution for the LC Green function is known only for the 
simplest form of the dipole cross section cx^ty) oc r\. With a realistic form of ^(ty) it 
is possible only to solve this problem numerically, what is still a challenge. Here we use 
the approximation suggested in ]18[ to evaluate the corrections arising from the finiteness 
of l c = 2vjM\ by multiplying the cross sections for incoherent and coherent production 
evaluated for l c — > oo by a kind of formfactor F mc and F coh respectively: 



a 



a 



[s,Q 2 ) =► a^ x (s,Q 2 )-F mc (s,l c (s,Q 2 
[s,Q 2 ) =► a^ A (s,Q 2 )-F c ° h (s,l c (s,Q 2 



where 



F mc (s,l c ) 



d 2 b dzp A (b,z) \F 1 (s,b,z)-F 2 (s,b,z,l c )\ 2 /(...)\ lc 



jcoh , 



S,lc) 



d 2 b 



-oo 
oo 



FAs, b,z) 



F 2 {s, b,z,l c ) 



dzp A {b,z)F l {s,b, z) e iz/l 

i i 7 

exp I --a 9N (s) J dz' p A (b,z' 

^ z 

-ct*jv(s) / dz' p A (b, z')Fi{b,z')e 



/(•••)l 



L=oo 



i(z—z')/l c 



(43) 
(44) 



(45) 
(46) 
(47) 
(48) 
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For the charmonium nucleoli total cross section cr^jv(s) we use our previous results [T|], 
essentially the expectation value of the dipole cross section cr gg -(r T , s) taken between the 
LC wave functions for the charmonia. The "formfactors" F inc for incoherent J/ip and ip' 
production are depicted in Fig. RL The observed nontrivial energy dependence is easy to 
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Figure 8: l c corrections for incoherent production of J/ip and if)' on C (thin 
curves) and Pb (thick curves) for the GBW parameterization of the dipole cross 
section a q g. 

interpret. At low energies l c <C Ra and the photon propagates without any attenuation 
inside the nucleus where it develops for a short time t c a cc fluctuation which momentarily 
interacts to get on mass shell. The produced cc pair attenuates along the path whose length 
is a half of the nuclear thickness on the average. On the other hand, at high energies l c 3> Ra 
and the cc fluctuation is developed long before its interaction with the nucleus. As a result, 
it propagates through the whole nucleus and the mean path length is twice as long as at low 
energies. This is why the nuclear transparency drops when going from the regimes of short 
to long l c . 

One can simplify expression (|45D assuming ct*jvTa(6) 1 which is a rather accurate 
approximation for J/i/j. Then one can expand the exponentials in (|47D-([48|) and obtain 0, 



Here 



F 2 A {<1) 



1 



l + F\{q c ) a mN {s)(T A ) 



1 - CF J/i)N (s){T A ) 



(Ta) 



d 2 b 



dz pA{b, z) e 



iqz 



(49) 



(50) 



is the so called longitudinal nuclear formfactor. 
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At low energies, q c = l/l c ^> 1/Ra, the formfactor F\(q c ) = and both the numerator 
and denominator in ( ^9|) decrease with energy due to the rising ctj/^n(s), but the denominator 
does it faster. This is why the ratio Eq. ([49]) rises with energy unless it drops due to the rise 
of the formfactor F^(q c ) ^ 0. The transition energy regime is related to the variation of the 
formfactor between and 1, i.e. at 



m N { 



{M 2 JM> + Q 2 )Rf/V3 = -so 



10 



(51) 



This energy increases with Q 2 and nuclear radius in accordance with Fig. 
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Figure 9: l c corrections for coherent production of J/if> and if)' on C (thin curves) 
and Pb (thick curves) for the GBW parameterization of the dipole cross section. 
Contrary to the case of incoherent production results for J/ip and if)' are very 
similar. 

Note that the Glauber approximation used to calculate F mc is rather good for J/tjj, but not 
for if)'. One could improve it applying the rigorous description based on the light-cone Green 
function formalism developed in || for exclusive production of vector mesons. However, 
it was done only for the unrealistic form of the dipole cross sections a q g cx r\ and for an 
oscillatory wave function of the vector meson. It is still a challenge to make this approach 
realistic. 

The results for the energy dependence of the formfactor F coh for coherent scattering are 
shown on Fig. |9|. In the limit of low energies l c — > the strongly oscillating exponential phase 
factor in ( fj6|) makes the integral very small and thus FJj^ w 0. Then the cross section rises 
with l c unless it saturates at l c ^> Ra when the phase factor becomes constant. Apparently, 
this transition region is shifted to higher energies for larger Q 2 and nuclear radius as is 
confirmed by the curves in Fig. || 
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5 Final results and conclusions 



Combining the results of the previous sections (i.e. including finite coherence length and 
gluon shadowing) we obtain the final results for ratios the Rj^ and R^r. Figs. |10| and [TT] 
show the s and Q 2 dependences for incoherent and Figs. and [13] for coherent production 
of charmonia. 
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Figure 10: Ratios R^ c for J/ip and ip' incoherent production on carbon, copper 
and lead as function of and at Q 2 = calculated with GBW parameterization 
of a q q. Solid curves display the modifications due to the gluon shadowing and 
finite coherence length l c while the dashed lines are without (same as on Fig. [I]). 

We summarize, in this paper we provide predictions for nuclear effects in exclusive elec- 
troproduction of charmonia on nuclei with realistic light-cone wave functions of charmonia 
and with dipole cross sections which have been already tested in a calculation of elastic 
charmonium electroproduction off protons. Such calculations still can be performed only at 
rather high energies where the coherence length exceeds the nuclear size (l c 3> Ra)- At these 
energies gluon shadowing which is absent at lower energies (l c ~ Ra) becomes an impor- 
tant and sometimes the dominant effect. Since no reliable experimental information about 
gluon shadowing is available so far, we have performed calculations of this effect employing 
the same light-cone dipole approach. We found sizable, but not large corrections (10-20%) 
which, however, keep rising with energy. It turns out that the magnitude of gluon shadowing 
is suppressed by a nonperturbative interaction of the light-cone gluons which is well fixed 
by data on large mass diffraction. The effects related to finiteness of the coherence length 
at somewhat lower energies were estimated too. 

We have calculated both coherent and incoherent cross sections as function of energy, 
Q 2 and momentum transfer. Our predictions can be tested in future experiments at high 
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Figure 11: Ratios R l £ c for J/ip and ip' incoherent production on carbon, copper 
and lead as function of Q 2 + M\ at fixed s = 4000 GeV 2 for GBW parameteri- 
zation of the dipole cross section a qq -. Solid curves correspond to the final result 
(with gluon shadowing and l c corrections) while thin curves are without (same 
as on Fig. 

energies with electron-nuclear colliders (eRHIC) . It also can be produced in peripheral heavy 
ion collisions (RHIC, LHC), and our predictions will be published elsewhere. 
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Appendix A Resummation expressions 

Matrix elements of the dipole cross section <j qq for the charmonia production including the 
effects of spin rotation were calculated in |l| for the specific form of the dependence 
provided by the KST and GBW parameterizations: 



L \ ' o / J 

In the nuclear case one needs to calculate matrix elements of the operators with a qq in 
the exponent. For incoherent and coherent scattering one can use Taylor expansion, for 
example: 




(A.l) 




(A.2) 
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1— f + 6 



+ 



(A.3) 



where f n =< cr n+1 > / < a > and c doesn't depend on initial and final states. So matrix 
elements of powers a q q n have to be calculated. The form ( |A.1|) allows to express any power 
of u q q{rT) as a sum of exponentials exp (— r\jr§) with reduced r\ =>- r 2 /m, i.e. <cr n > can 
be expressed in terms of <a> with reduced r^. 

The rate of convergence of such series slows down with increasing a (larger s) and increas- 
ing nuclear thickness (c oc T^(6)). In our calculations we used a more efficient expansion, 
which contains matrix elements < o n > in the exponent. This form can be obtained when 
applying the identity 



l + x = e log{1+x) = exp (x - y + y + 



(A.4) 



cxe 



to ( |A.2| ) (with x = cfi + ■ • •) and ( A.3| ) (with a; = — c/i/2 + ••■). As a result one can get 

= < <r > • exp (-ch + j (h - ft) + ■ ■ -J , (A.5) 



1 -e~ ca 

which converges noticeable faster. Fig. [14| shows that while one needs about six terms (n = 6) 



to obtain an accurate result using Taylor expansion ( |A.2|) , one has already satisfactory results 
for n = 2 in the method, where one resums in the exponent. 
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